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Abstract. This is intended to be a brief introduction and overview of Hofava-Lifshitz gravity. 
The motivation and all of the various versions of the theory (to date) are presented. The 
dynamics of the theory are discussed in some detail, with a focus on low energy viability and 
consistency, as these have been the issues that attracted most of the attention in the literature 
so far. Other properties of the theory and developments within its framework are also covered, 
such as: its relation to Einstein-aether theory, cosmology, and future perspectives. 



1. Introduction 

1.1. Gravity and renormalizability 

It is a well known fact that general relativity is not a renormalizable theory. This means that, 
successful as it may be as a classical theory of gravity, it should be viewed as an effective theory 
that breaks down at some scale. Beyond that scale general relativity is not enough to describe 
the gravitational interaction or spacetime itself (depending on the perspective) and one cannot 
construct its quantum counterpart using conventional quantization techniques. 

On the other hand, given that within this perspective general relativity is an effective theory, 
the Einstein-Hilbert action contains only the lowest order terms in a curvature expansion. A 
reasonable question is then whether a renormalizable theory can be obtained by including higher 
order curvature terms. Intuitively this could work, as such terms would modify the propagator of 
the graviton at high energies. Indeed, in 1977 Stelle showed that theories whose action includes 
invariants quadratic in the curvature are renormalizable [1]. However, unfortunately this comes 
at a very high price, as such theories contain ghost degrees of freedom and are, therefore, not 
unitary. This is because by adding higher order curvature invariants in the action, one adds 
higher order time derivatives. 

One could attempt to modify the propagator by adding higher order spatial derivatives 
without adding higher order time derivatives. This could presumably lead to a theory with 
improved ultraviolet (UV) behaviour without having to face the dreadful consequences of having 
higher order time derivatives. Obviously, following such a prescription requires treating spatial 
and time derivatives on a different footing and as such is in clash with Lorentz invariance. On 
the other hand, since the modified behaviour of the propagator is strictly needed in the UV, it 
is conceivable that Lorentz invariance could be recovered in the infrared (IR), or at least that 
Lorentz violations in the IR could stay below current experimental constraints. 

Hofava's proposal [2, 3], which has led to a spree of publications recently and is now commonly 
referred to as Hofava-Lifshitz gravity, can be seen as an attempt to put these heuristic arguments 
in a rigorous framework where one can test whether they are indeed valid when it comes to 



gravity theories. In what comes next, an overview of this framework is given and its major 
developments to date are briefly discussed. This is far from an exhaustive account of an already 
very extensive literature, but hopefully it covers the most critical topics and provides a concise 
introduction to the subject. 1 



1.2. Lorentz violations as afield theory regulator 

What has been discussed above is essentially the idea of giving up Lorentz invariance in order 
to keep a quantum field theory finite. Clearly, this idea goes beyond gravity theories and it has 
been considered in the past for other fields, see for example Ref. [5]. In fact, it would be best to 
start from a simpler example, such as a scalar field. This section will closely follow Ref. [6]. 

Let us first of all point out that there is nothing wrong with using Lorentz violations as 
a field theory regulator from a perspective of logical consistency, as there is no reason why a 
theory should necessarily exhibit Lorentz symmetry in the far UV. On the other hand, Lorentz 
violations are severely constrained in a wide range of energies and especially in the IR. Therefore, 
as mentioned before, the real question is whether one can construct a field theory that exhibits 
Lorentz violations which in the far UV lead to renormalizability, but remain suppressed below 
experimental accuracy at lower energies. This is difficult to answer without concrete models and 
enough to motivate further investigation. 

Consider the scalar field action 
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= J dtdx d \ <P 2 - a m 0(- A) m + ]T g n ^ I , (1) 
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where ' = d/dt, A = V 2 is the spatial Laplacian and z and M are positive integers still to be 
specified. A theory is said to be "power counting renormalizable" when all of its interaction terms 
scale like energy to some non-positive power, as in this case Feynman diagrams are expected 
to be convergent (or have at most a logarithmic divergence). To check whether this is true 
or not we first have to choose the engineering dimensions of space and time. The engineering 
dimensions 

[dt] = [ K ]' z , [dx] = [k]-\ (2) 

where k is a place holder symbol with dimensions of momentum, are compatible with the choice 
a z = 1 and the requirement that the action has to be dimensionless. This choice of engineering 
dimensions, therefore, appears to be natural in the UV. The dimensions for the scalar field then 
turn out to be 

[4>] = [4 d - z)/2 - (3) 

It remains to check the scaling for the rest of the terms. Requiring that an interaction term scales 
like energy to some non-positive power is equivalent to the requirement that the coupling of this 
interaction scales like energy to some non-negative power (since the action is dimensionless). 
Also, we find it convenient to work in terms of momenta instead of energies. So, our theory will 
be power counting renormalizable if all couplings have non-negative momentum dimensions. 
It is straightforward to verify that 

M = M 2(z_m) , [9n] = [«]^-»("-*)/2. ( 4 ) 

It is then obvious that a m has non negative momentum dimension for all values of m. It is also 
easy to realize that, for z > d, g n has non- negative momentum dimensions for all values of n. 
When z < d, g n has non-negative momentum dimension only when n < [2(d + z)]/(d — z). One 



1 See also Ref. [4] for an early review of the projectable Hofava-Lifshitz gravity. 



can then easily recover well-known results. For example one can verify that the case d = 3, 
z = 1, M = 4 which corresponds to the usual relativistic <p A theory in 3 + 1 dimensions is a 
power counting renormalizable theory. The d = z = 1 case, i.e. the relativistic scalar in 1 + 1 
dimensions, is power counting renormalizable for all values of M. In fact this case belongs to 
the class z = d and is directly analogous to the case d = z = 3, which will also have the same 
property and is of interest to us. In fact one can argue that any theory with z = d is actually 
finite with normal ordering, while theories with z > d are finite even without normal ordering, 
see Ref. [6] for a more detailed discussion using the "superficial degree of divergence" . 

How would things go if instead of a scalar field we wanted to consider a graviton? The main 
difference lies in the graviton self-interaction vertices. In the case of the scalar field considered 
above, momenta did not enter these interactions. For a graviton, however, one would have to 
deal with self-interaction vertices which contain spatial derivatives [7]. This introduces some 
further complication, but it is not enough to spoil the power-counting renormalizability, as long 
as z > d, i.e. as long as the action contains operators with at least 2d spatial derivatives. The 
presence of self-interaction vertices which contain spatial derivatives is enough though to argue 
that, unlike the scalar theory consider above, its graviton analogue will not be finite (even with 
normal ordering) [7]. 

All of the arguments presented here support the statement that field theories that contain 
at least 2d spatial derivatives in d + 1 dimensions are power counting renormalizable. Power 
counting renormalizability is a strong indication that a theory is indeed renormalizable. Little 
progress has been made so far in the case of gravity to check renormalizability beyond power 
counting for such theories. 

2. From the Lifshitz scalar to non-relativistic gravity 

We proceed to consider the exact structure of a gravity theory with the characteristics mentioned 
in the previous section. We will work in 3 dimensions from now on but generalizing to more 
than 3 dimensions is straightforward. The basic requirement is that the theory should have only 
two time derivatives but at least 2d spatial derivatives, that is 6 in our case. The fact that we 
want to have more spatial than time derivatives implies that we need to treat space and time on 
different footing. This is naturally achieved by working in the Arnowitt-Descr-Misncr (ADM) 
decomposition of spacetime 

ds 2 = -N 2 c 2 dt 2 + g ij (dx i + N { dt)(dx j + N j dt). (5) 

One is essentially forced to pick a preferred foliation of spacetime in order to write down the 
action according to the prescription laid out in the previous section. Consequently, such an action 
cannot be invariant under standard diffeomorphisms, as general relativity is. It can, however, 
be invariant under a more restricted set: foliation preserving diffeomorphisms, i.e. space- 
independent time reparametrization together with time-dependent spatial diffeomorphisms 

t->i(t), x i ^x i {t,x i ). (6) 

This is the symmetry that the action will have to respect. 

It is straightforward to realize that the only covariant quantity under spatial diffeomorphisms 
that contains a time derivative of the spatial metric is the extrinsic curvature 

K H = ^ {9ij - ^ - VjNi} , (7) 

which also transforms like a scalar under time reparametrization. A dot denotes differentiation 
with respect to the time coordinate and Vj is the covariant derivative associated with the spatial 



metric gij. We want the theory to be second order in time derivatives, so we need to consider 
terms quadratic in the extrinsic curvature. Note that there are no invariants under the symmetry 
mentioned above that one can construct with time derivatives of the lapse N or the shift Nj, 
without including higher order time derivatives of as well. Therefore, the most general action 
one can write is 

M 2 r 

S =^rJ d 3 xdtN^g{K^K i:j - XK 2 - V( 9ij ,N)} , (8) 

where M p \ is a constant which, with some foresight, we identify with the Planck mass, g is the 
determinant of the spatial metric g^ , and A is a dimensionless running coupling. V generically 
can depend on g$j and N and their spatial derivatives. It does not contain time derivatives, 
neither does it depend on the shift iVj, as the symmetry of the theory does not allow one 
to construct any suitable invariants. Finally, power counting renormalizability requires V to 
contain terms which are at least sixth order in spatial derivatives. From now on we will restrict 
ourselves to theories that contain sixth order time derivatives but not higher, as these are the 
simplest ones which satisfy the power counting renormalizability requirement. 



3. Potential, constraints and the various versions of the theory 

Clearly there are numerous terms one could include in V. Different choices have led to different 
versions of the theory. In this section we list these versions, briefly presenting their basic 
characteristics and differences. We refrain from discussing or commenting upon their viability 
or consistency, as this, together with their phenomenology, will be the subject of the coming 
section. 

3.1. Detailed balance 

Hofava proposed a symmetry that V would have to satisfy, that would drastically reduce the 
number of invariants one should consider [3]. This symmetry is dubbed detailed balance and 
it is inspired by condensed matter systems. It sums up to the requirement that V should be 
derivable from a superpotential W: 

V = E^g ijkl E kl , (9) 

W ^ ere 1 sw 

E T/gSgij 1 ^ 

and 

gijM = l ^ gikg jl + gilgjk^ _ Xg i jg M 

The most general action one can write with V satisfying the conditions above is 

M 2 , f { A- 4 2rv 2 ft f^ k ft 2 

S db = ^j d^xdtN^gl I&Kv - XK 2 - j^C^ + -^f—^VjR 1 , - j^R ij Rij 

/3 2 1-4A 2 /3 2 C 3/3 2 C 2 2 1 

+ ^T^x R+ Hx R -T=k M ^^ ^ 

where e l ^ k is the Levi-Civita symbol, 

C ij = ^=V* (V, - -R5\} (13) 



is the Cotton tensor, which in 3 dimensions plays the role of the Weyl tensor, and a, (3 and £ 
are dimensionless couplings. Notice that there are only 3 new couplings for a total of 6 terms 
in V. 

The advantages of detailed balance is that it drastically reduces the number of terms one needs 
to consider and that it introduces a superpotential, which might simplify quantization. On the 
other hand, the action above contains a term that violates parity (the fifth order operator). 2 
Additionally, as is obvious, the last term in the action, which plays the role of a cosmological 
constant, is restricted to have the wrong sign with respect to the one suggested by observations 
when A > 1/3, as pointed out in Refs. [9, 10] (A is supposed to be close to unity, if severe Lorentz 
violation are not to occur). 3 Note that there is nothing fundamental about detailed balance, so 
it is safe to consider it just a simplicity assumption. 

3.2. Projectable Hofava-Lifshitz gravity 

Another possible simplifying restriction proposed by Hofava in Ref. [3] is that of projectability, 
which is the assumption that the lapse is just a function of time. That is A = N(t). Again there 
is no fundamental principle behind such an assumption. Hofava's main motivation for initially 
considering it was that under this assumption one has enough gauge freedom to set N = 1, as 
in general relativity. This is not possible without projectability, as the symmetries of the action 
allow only space independent time reparametrizations, as explained earlier. 

In any case, projectability drastically reduces the number of invariants V can include even 
if detailed balance is not imposed, as all terms with spatial derivatives of N would now vanish. 
Let us see that in more detail, following the line of Refs. [9, 10]. Since N is now just a function 
of time there is no invariant one can built out of it without using time derivatives. So V will 
only depend on the metric and its spatial derivatives. Essentially, this means that the action 
should include all of the curvature invariants one can construct with gij with up to six spatial 
derivatives. 

An important simplification comes from the fact that we are in 3 spatial dimensions. 
Therefore, the Weyl tensor vanishes identically and the Riemann tensor can be expressed in 
terms of the Ricci tensor. Taking Bianchi identities into account as well and ignoring surface 
terms, one can arrive at the conclusion that the most general action is 

M 2 r [ 

S p = -f- J SxdtN^l K^Kij - XK 2 - g M 2 , - 9l R - g 2 M~ 2 R 2 - g 3 M~ 2 

-g 4 MpJ 4 R 3 - g 5 M pl A R(RijR ij ) - g 6 R l jW k R k t 

-g 7 M p - 4 RV 2 R - g 8 M pl 4 ViR jk V l i?^| , (14) 

where the gi are dimensionless couplings. As long as we have not coupled matter to the theory 
we are free to rescale the coordinate in order to set g\ = —1, which is the value it has in 
general relativity. We then have 8 remaining couplings gi which can be used to tune the scales 
suppressing the various operators. The following remarks are in order: 

• We have suppressed parity violating terms. 

• Enforcing projectability on top of detailed balance (without suppressing parity violating 
terms) would leave action (12) unaffected, apart from the fact that N would become just a 
function of time. As we will see later, this is quite an important difference. 

2 Parity violating terms, within or outside the framework of detailed balance, have been shown to lead to chiral 
primordial gravitational waves [8]. 

3 This is a bare cosmological constant. Also, one could analytically continue the parameters as explained in 
Ref. [11], but in this case W would cease to be real. 



• There are just 3 more operators in the most general projectable action than in the one with 
detailed balance. In this sense, detailed balance does not bring significant simplification 
once project ability has been assumed. For this reason, and for the drawbacks listed in the 
end of the previous section, we will not consider detailed balance further when assuming 
projectabiltity. 

• go controls the value of the cosmological constant, which, unlike the case with detailed 
balance, is not restricted. 

• There are two types of Lorentz violating terms in the action. The ones in V, which are 
suppressed by some scale that can be determined by tuning the couplings §2 to g% , and one 
that comes from the kinetic part due to the fact that A is not necessarily equal to 1. This 
will be a generic feature of all versions of the theory as we will see. Clearly, the term of the 
second kind is far more dangerous, as it introduces Lorentz violations at low energy scales. 

• All couplings are running. One can, therefore, hope that A will run to 1 in the IR (or 
sufficiently close to it for experimental constraints to be satisfied), the rest of the operators 
will be heavily suppressed, and the theory will effectively reduce to general relativity. 
Diffeomorphism invariance and Lorentz invariance would emerge as IR (approximate) 
symmetries. Whether or not this will indeed be the case requires a study of the 
renormalization group flow which is still pending. 



3.3. Non-projectable Hofava-Lifshitz gravity 

We now proceed to consider the case where neither detailed balance nor projectability are 
enforced. This version of the theory is called non-projectable. To avoid confusion it is worth 
pointing out explicitly, that the version with detailed balance can be a non-projectable version. 
However, in this section we wish to go beyond that. Once one abandons detailed balance, 
adding just some specific choice of extra terms is not really an option. Radiative corrections 
will generate all possible terms compatible with the symmetries of the theory and, thus, all such 
terms should be taken into account, much like the projectable case treated above. However, as 
first pointed out in Ref. [12], it is not only curvature invariants of gij that one can include in V 
in this case. One can also use the quantity 

ai = dilnN, (15) 

as contractions of it with itself or curvature terms also lead to invariants. The lowest order 
invariant one can construct with ctj is a l ai, which comes at the same order as R. The action will 
then be of the form 

M 2 r ( 1 1 ) 

S np = J d 3 xdtN^ I K^Kij -\K 2 + CR + r] a t a l + ^ L 4 + \ , (16) 

where L4 and Lq include all possible 4th and 6th order operators respectively that one can 
construct using aj and g^. For example, R 2 , (a J aj) 2 , V*Vj-R and a^ajR 1 ^ are some of the terms 
that are included in L4. 4 As was the case before for gi, £ can now be set to unit by a coordinate 
rescaling when matter is not coupled to the theory. This version of the theory was first studied 
in Ref. [12] in order to resolve some inconsistencies in the dynamics of the non-projectable 
version with detailed balance, as we will discuss shortly. Although it is often referred to as an 
"extension" of Hofava-Lifshitz gravity, this terminology is probably unfortunate: this is Hofava- 
Lifshitz gravity, i.e. a theory constructed consistently according to the prescription described in 
Hofava's paper Ref. [3] and without any restriction, such as detailed balance or projectability. 

The following points are worth stressing for the non-projectable version (apart from some of 
the remarks already made for the projectable case and apply here as well). 



4 V l Vi7? and other similar terms are not surface terms, as the N in the 4-volume is space dependent. 



• For simplicity the cosmological constant term has been suppressed, but it can be 
straightforwardly restored. 

• The scales Ma and Mb suppressing the higher order operators have been left arbitrary. 
This is in analogy with having arbitrary dimensionless couplings gi in the projectable case. 

• The number of operators in the non-projectable case is an order of magnitude larger than 
that in the projectable case. 

• If general relativity were to be recovered in the IR, not only A would have to run to 1 but 
also 7] would have to run to 0. 

3.4- Imposing further symmetries 

All of the actions presented above, which correspond to different versions of Hofava-Lifshitz 
gravity, are invariant under foliation preserving diffeomorphisms only, and not the full set of 
diffeomorphisms. So, even though the action is in all cases quadratic in the time derivatives of 
gij, less symmetry generically implies that more degrees of freedom will be excited. We will see 
explicitly in the next section that this is indeed the case, and that this can, in most cases, have 
undesirable consequences for the consistency and viability of the theory. 

A way to do away with the extra degrees of freedom is to attempt to nontrivially extend the 
gauge symmetry of the theory, so that it will have as many generators per spacetime point as 
general relativity. It was pointed out in Ref. [13] that this could be done if the action would be 
suitably modified to acquire an extra local U(l) symmetry. More precisely in Ref. [13] a step 
by step construction of the action led to 

M 2 r 

S extraU(1) = J dPxdtNy/g^Kij — K 2 — V{ 9ij ,N) 

+v® ij {2Kij + ViVji/) - A(R - 2Q)/n} , (17) 
where Q is a constant, A acts as a Lagrange multiplier, v is an auxiliary scalar field, 

QV = R ij -IgVR + ngV, (18) 

and V should as usual contain 6th order operators and may or may not satisfy detailed 
balance. Note that A transforms as a spatial scalar and a time vector under foliation preserving 
diffeomorpisms, that is 

A^A + fA + fA + ediA, (19) 

where is the generator. On the other hand, under a local U(l) gauge transformation with 
generator if, A and v transform as 

A -»• A + ip - N*Vi(p, (20) 
+ <p- (21) 

See Ref. [13] for more details on constructing the action, as well as for discussions on 
anisotropic scaling, hamiltonian formulation, linearization, etc. It was argued there that in 
the IR this theory has very good chances to reduce to general relativity, as the high order 
operators in V would be suppressed (as in the previous versions of the theory), the extra fields 
in the action are just auxiliary fields, and the extra gauge symmetry leads to just the usual 
spin 2 graviton excitation. Moreover, it was claimed in Ref. [13] that the U(l) symmetry forces 
the coefficient A of K 2 in the action, which was a running coupling in all other versions of the 
theory, to be equal to 1 here, as in general relativity. 



This last claim has been contested recently in Ref. [14] where it has been claimed that the 
extra U(l) symmetry does not necessarily require A = 1 and that even for arbitrary A the theory 
propagates only a spin 2 graviton. If this is true, A has to remain a running coupling and flow 
close to 1 in the IR for general relativistic phenomenology to be recovered according to Ref. [14]. 
So, at the moment, whether or not the value of A is fixed to 1 can be considered a matter of 
debate. 

Another potential difficulty with this version of Horava-Lifshitz gravity has to do with matter 
coupling. It is straightforward to notice that the Lagrange multiplier in action (17) forces R to 
be constant. Even though this is not necessarily a problem in vacuum, it cannot be generically 
the case once matter is coupled to the theory. This implies that matter would have to be coupled 
to A if the theory is to resemble general relativity at low energies, and that a suitable coupling 
that leads to standard phenomenology in the IR should exist. Whether this is indeed the case 
requires further investigation (some first steps where done in Ref. [14]). It is also worth pointing 
out that action (17) is susceptible to radiative corrections which could modify its form. 

Since this last version of Horava-Lifshitz gravity has been proposed very recently and has 
not yet been studied extensively, we will not consider it further. Its cosmological aspects where 
recently studied in Ref. [15]. 

4. Projectable version: Dynamics, consistency and low energy behaviour 

4-1. Degrees of freedom and propagators 

We return now to the projectable version of the theory and take a closer look at its dynamics. 
We will consider the most general action in this version, action (14), i.e. we will not enforce 
detailed balance but our results will include it as a subcase. The field equations for this action 
have been derived in detail in Ref. [10]. We will refrain from rewriting them here due to space 
limitations but we will make the following remark: since N = N(t), the variation with respect 
to N will not lead to a usual local super-Hamiltonian constraint, but to a global Hamiltonian 
constraint. That is, the constraint will come in the form of an integral over space. This will be 
a major difference with respect to the non-projectable case. 

We now move straight to the linearization. If for simplicity one sets the cosmological constant 
term to 0, then flat space is a suitable background. After suitable gauge choices, see Ref. [10], 
one finds that the theory propagates a spin-2 mode, the usual graviton, that now satisfies a 
modified dispersion relation, i.e. 



where Hij is transverse and traceless. Notice that the couplings g% and gs control the scale 
that the lorentz violating terms become important in the dispersion relation. Recall also that 
in vacuum one has the freedom to rescale the coordinates and set g\ = — 1. 

As mentioned earlier, less symmetry generically means more degrees of freedom, so since 
the theory is invariant under foliation preserving diffeomorphisms only, we expect to find extra 
excitations. This is indeed the case. There is an extra scalar degree of freedom, whose linearized 
dynamics are governed by the action [10] 



= - gid 2 + g 3 M- 2 d 4 + g s M^dP H ijt 



(22) 
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(23) 



where 




1-A 
3A - 1 



(24) 



and we have set g\ = — 1 by rescaling the coordinates, in order to guarantee the stability of the 
spin 2 graviton. Given the overall minus sign in eq. (23), the scalar mode is a ghost whenever 
1 > A > 1/3. On the other hand, the scalar is classically unstable at low energies whenever 
A>lorA<l/3 [10]. It is conceivable that the instability can be cut off by the higher order 
derivatives. If we call M* the scale where the higher order term would take over, then the time 
scale of the instability would be at least l/\ch\M*. The relation between M* and M p \ is clearly 
governed by the g. t coefficients in eq. (23). On the other hand, if the instability is not to develop 
within the lifetime of the universe then \ch\M* < Hq, where Hq is the Hubble parameter. [16]. 
Experiments rule out modifications of Newton's law at scales above Wfim which corresponds 
roughly to the constraint M* > O.leV. In turn, this implies |1 — A| < 10~ 61 . This value is clearly 
very low and it seems unlikely that the renormalization group flow could drive A to it. 

In Rcfs. [17, 18] linearization around a de Sitter, instead of Minkowski, background was 
considered. As expected, the large scale modes exhibit better behaviour. However, qualitatively 
the result is the same overall: A needs to be sufficiently close to 1 for the instability not to be 
a concern. Note that requiring A to be very close to one is not only a potential fine tuning 
problem. Ch — > as A — > 1, which means that the low momentum phase velocity of the scalar 
mode will be much smaller than the speed of light for such values of A. When a velocity of a 
certain mode is smaller than the speed of light in a medium, relativistic particles traveling in 
this medium decay via the Cherenkov process [19]. The fact that we observe high energy cosmic 
rays which need to travel a significant distance to reach us is evidence enough that such modes 
do not exist. 



4-2. Strong coupling 

In the previous section we only considered the linearized dynamics as described by the part of 
the perturbative action quadratic in h. However, it has been pointed out by several authors 
that such a perturbative treatment breaks down when A approaches 1 and the scalar mode gets 
strongly coupled [20, 21, 16]. To see how this comes about we give below the cubic interactions 
of h, 



M, 



pi 



dtofx \ h(dh) 2 - ^rhdih^h + - 
ci d 2 2 



did. 



1 h 



3 



h) - 



(H + 1) 



hh 2 



(25) 



where we have neglected cubic interactions coming from higher order operators. The redefinitions 

—1/2 

t = \Ch\t and h = c h M p \h canonically normalize the lower operator part of the quadratic 
action (23) (first two terms). After these redefinitions the cubic action reads 



M 3 / 2 M pl 



/ 



did 3 x { c\ h{dhf 



2h'dih-,h' 
d 2 



d 2 n 



- (2c 2 + l)h{h>) 



(26) 

where now ' = d/df. As is obvious from the equation, there are cubic interactions which are 
suppressed by the scale Ic/J^Mpi with respect to the quadratic ones. Therefore, the theory 
gets strongly coupled at the scale M sc = \ch\ 3 ^ 2 M p \. 

Given the constraints on A from stability discussed previously, the strong coupling scale is 
phenomenologically unacceptably low, as we know that we can treat gravity perturbatively at 
low energies. There is also the issue of renormalizability: our arguments for power-counting 
renormalizability given above where based on the validity of the perturbative treatment at all 
energies. If there is strong coupling such arguments simply fail. 

The only subtlety here is the fact that the strong coupling arguments presented so far neglect 
the higher order operators. However, for such operators to become important and invalidate 



our treatment, the strong coupling scale M sc would have to be higher than the scale M* that 
suppresses them. Given how low M sc is here, this would imply that higher order operators 
modify the graviton dynamics at very low energies, which would be in conflict with current 
experimental evidence. 

Since the analysis presented here is linear, one could ask whether including non-linear effects 
could resolve the strong coupling problem, by an analogue of the Vainstein mechanism in massive 
gravity [22], i.e. a non-perturbative restoration of the A — > 1 limit. Recently, in Ref. [23] it has 
been claimed that this is indeed the case in spherically symmetric, static configurations. This 
matter deserves further attention. 



5. Non-projectable version: Dynamics, consistency and low energy behaviour 

5. 1 . Degrees of freedom and propagators 

We now turn our attention to the non-projectable version of the theory. As before, we will 
not enforce detailed balance in order to obtain the most general results possible. However, we 
will discuss explicitly that case as well. The dynamics of the spin 2 graviton are essentially the 
same as in the projectable version so we will not consider them again in more detail. As in the 
projectable case, also here there is an extra scalar degree of freedom, and we will focus on that. 
The low energy dynamics of this scalar graviton are governed by the action 



S^> = -M p 2 ! J d 3 xdt 



lh 2 + ^hd 2 h 
c h v 



(27) 



which one obtains by linearizing action (16) around flat space and considering only the lowest 
order operators. £ has been set to 1 by a suitable coordinate rescaling. Note that Ch is defined 
in the same way as above and, thus, it is not the low momentum phase velocity of the scalar in 
the non-projectable version. The square of the latter is now 

c' h 2 = cl V -^. (28) 

For h not to be a ghost one needs c\ < 1, and for h to be classically stable one needs c'j* > 0. 
These conditions are satisfied whenever [12] 

A > 1, < n < 2, (29) 



or 

A < 1/3, < t) < 2. (30) 

In the second region in the parameter space A is far from 1, which is the value it has in general 
relativity, so we will not consider this option further. The presence of the operator a^a* in action 
(16), which is a lowest order operator and contributes to the quadratic action, has drastically 
improved the behaviour of the scalar graviton with respect to the projectable case we examined 
previously. In particular, there is now a region in parameter space, described by the constraints 
(29), for which h is neither a ghost nor classically unstable [12]. Even though we have neglected 
higher order operators here, the behaviour we have found would remain qualitatively the same 
if we had included them. 

Let us examine what would happen if we had enforced detailed balance. The a l ai term in 
the action would not be allowed which means r/ = 0. In this case, the coefficient of the second 
term in (27) blows up and the scalar field appears to freeze. Things, however, are actually more 
complicated, as a more detailed analysis shows [21]. When detailed balance is enforced, the 
scalar does propagate around backgrounds that are neither static nor homogeneous. In fact it 
appears to satisfy a first order differential equation, which constitutes a serious problem for the 



theory, see Ref. [21] for more details. Similar conclusions have been obtained by considering 
the hamiltonian formulation of the theory with detailed balance [24, 25]. In particular, the 
Hamiltonian constraint is not automatically preserved by time evolution (second class constraint) 
and the theory propagate 5/2 half degrees of freedom, the 1/2 corresponding to the scalar mode 
which satisfies a first order equation and the rest to the two graviton polarizations. Such a 
theory, though mathematically consistent, exhibits no time evolution and is, therefore, physically 
meaningless [25]. 



5.2. Strong coupling 

We just saw that, judging from the quadratic action, the non-projectable version exhibits 
improved scalar dynamics with respect to the projectable version, as long as detailed balance 
is not enforced. However, we found earlier that the scalar mode in the projectable version gets 
strongly coupled at unacceptably low energies, so one has to check whether this is the case in 
the non-projectable version as well. To do so, we write down the cubic action for the scalar [26]: 
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Same as before for the projectable case, we canonically normalize the low energy quadratic 
action. This requires the redefinitions 



V \Ch\ 
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(32) 



Recall that we are only interested in the region of the parameter space where conditions (29) 
are satisfied. Under the redefinitions the cubic action reads 
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where we have expressed everywhere Ch in terms of the physical quantity c' h . From eq. (33) one 
directly infers that the cubic interactions are suppressed with respect to the quadratic ones by 
various scales /(|A — l\,r])M p \, where / is an algebraic function whose functional form depends 
on which terms one considers. The scalar mode would become strongly coupled at the lowest of 
these scales. 

Both | A — 1 1 and r\ essentially measure deviations from Lorentz invariance which are present at 
arbitrarily low scales. They are to be small for the theory to avoid experimental constraints on 
Lorentz violations. It is rather straightforward to verify that in this case, also the scale that the 
theory gets strongly coupled has to be low. Instead of performing a general analysis leaving both 
| A — 1 1 and n small but arbitrary, we will focus on the most interesting case where d h ~ 1. This 
value of c' h is the one dictated by Cherenkov radiation constraints [19], see discussion above. We 
then have n ~ | A — 1| and the strong coupling scale turns out to be M sc ~ ^/fjM p \ ~ y^A — l|M p i 
[26, 27]. 

Since rj, |A — 1| <C 1 the strong coupling scale will be much lower than the Planck scale. 
Absence of preferred frame effects in the Solar system observations requires rj, |A — 1| < 10~ 7 



[28], which in turn implies that M sc < 10 15 GeV. Therefore, the strong coupling scale is too 
high to be phenomenologically accessible from gravitational experiments. In this sense, strong 
coupling in non-projectable Hofava-Lifshitz gravity would pose no threat if the theory were to 
be treated as an effective theory describing gravity at low energies. On the other hand, the 
mere presence of strong coupling even at relatively high energies, as in this case, is enough to 
seriously question the power counting renormalizability of the theory [26]. Indeed, as mentioned 
also earlier when considering the projectable case, renormalizability arguments were based on 
the assumption that a perturbative treatment can be used to arbitrarily high energies. 

However, the strong coupling arguments given here neglect the role of the higher order 
operators in the action. As proposed in Ref. [29], a way to avoid strong coupling is to lower 
the scale that suppresses the higher order operators in action (16) below M sc , i.e. impose that 
Ma ~ Mb ~ and M* < M sc . This clearly requires the introduction of a second scale, 
different from M p \, as well as a peculiar hierarchy of scales. Also, it does seem to involve some 
large dimensionless couplings as M p \ 3> M*. However, as it has been argued in Refs. [29, 30] 
this hierarchy of scales is technically natural. 

An important consequence of this effective constraint on M+ coming from the requirement 
to avoid strong coupling, first pointed out in Ref. [26], is that now M* becomes bound from 
both above and below. Being the scale that suppresses higher order derivatives in the dispersion 
relations, one would expect that the larger M* is, the better. Assuming that this same scale 
suppresses higher order derivatives in the matter dispersion relations (which would also be 
modified in a way similar to gravity), and in particular photons, places a strong lower bound on 
M*: The absence of a time delay in different frequency 7-rays coming from distant astrophysical 
sources implies that [31, 32] M± > 10 11 GeV. However, now one has also an upper bound if 
strong coupling is not to be a problem, and the combined set of constraints yields 



There seems to be a comfortable window for M* within which non-projectable Hofava-Lifshitz 
gravity (without detailed balance or other restrictions) avoids strong coupling without exhibiting 
detectable Lorentz violations, at least with current experimental accuracy. Improving the 
accuracy with which we can constrain both preferred frame effects in the solar system and 
modifications in the dispersion relations could potentially close this window. 

5.3. Relation to Einstein- Aether theory 

Non-projectable Hofava-Lifshitz gravity is a theory with a preferred foliation, which will exhibit 
Lorentz invariance violations even at low energies. This should be clear by now by inspecting the 
action (16) and taking into account that, given the dynamics of the scalar discussed previously, 
one can no longer hope that A and 77 will run to 1 and respectively in the IR, in order for 
Lorentz symmetry to emerge. Note that the operators related to these couplings are of lowest 
order, i.e. dimension 2. 

A model theory for gravity with preferred frame effects is Einstein- aether theory [33, 34]. It 
is the most general theory for a unit timelike vector field coupled to gravity (but not to matter) , 
which is second order in derivatives. The most general action for Einstein-aether theory, up to 
total derivative terms and setting aside matter coupling, is 



where g^ u is the 4-dimensional metric, g its determinant, R the Ricci scalar of this metric, V 
is the covariant derivative associated with g^ u , 



10 15 GeV > M* > 10 n GeV. 



(34) 




(35) 



(36) 



with M a ?^ defined as 



M a ^ = c^sT + c 2 g^f u + c 3 g au 9^ + cfifvPg^. (37) 

Greek indices run from to 3, the c% are dimensionless coupling constants, and it is assumed that 
is constrained to be a unit vector, g^u^u^ = 1. This constraint can be explicitly imposed by 

the use of a lagrange multiplier, or implicitly taken into account in the variation of the action 

by allowing only variations that respect it. 

One may wonder whether there is any relation between Einstein-aether theory and the (naive) 

IR limit of non-projectable Hofava-Lifshitz gravity (that is, a theory described by action (16) 

without the higher order operators suppressed by Ma and Mb). It was shown in Ref. [35] that 

the latter is actually a limiting case of the former. Let us briefly go through the arguments 

presented there. 

Starting from the action (35), impose the restriction that the aether be hypersurface 
orthogonal: 

u a = — , 38 

y/g^TdvT V ; 

where the unit constraint on the aether has been taken into account. T will satisfy a dynamical 
equation. However, since it is a scalar, on shell its dynamical equation will be implied by the 
contracted Bianchi identity, provided that T is not a constant (this can be shown for a general 
scalar field coupled to gravity, provided that the theory is manifestly covariant and the other 
fields are assumed to satisfy their field equations). However, T cannot be a constant as it now 
defines the foliation. Thus, we are allowed to not explicitly impose the equation of motion for 
T, which in turn implies that T itself can be chosen as the time coordinate. We then have 



u, 



a = 5 aT (g TT r 1/2 = N5 aT . (39) 



This choice selects a preferred foliation with N being the lapse function of this foliation. 
Substituting eq. (39) into the Einstein-aether theory action (35), one gets after some algebra 
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(41) 



where we use the convention Cjj = q + cj. 5 Clearly, action (40) is the same as action (16) 
without the higher order operators. Therefore, the IR limit of non-projectable Hofava-Lifshitz 
gravity is equivalent to Einstein-aether theory with a hypersurface orthogonal aether. 

This equivalence provides new insights to both theories. For Einstein-aether theory it implies 
that there may exist a UV completion, at least for some subclasses. For non-projectable Hofava- 
Lifshitz gravity it tells us that it can be viewed as a theory were a preferred foliation is 
dynamically defined by a scalar field. Note that one would expect to be able to supplement 
Einstein-aether gravity (with a hypersurface orthogonal aether) with suitable higher order 
operators and construct a covariant equivalent of the full non-projectable Hofava-Lifshitz gravity 
(see Ref. [36] for an early attempt). 

Clearly, on account of the equivalence, some results obtained for Einstein-aether theory go 
through to the IR limit of Hofava-Lifshitz gravity and vice versa. A characteristic example 



5 Note that there is a typo in the formula relating £ and C13 in eq. (10) of Ref. [35]. 



is spherically symmetric solutions, as the aether is always hypersurface orthogonal in spherical 
symmetry [35, 30]. For example, some spherically symmetric solutions found in Ref. [37] for 
Hofava-Lifshitz gravity had actually already been found earlier for Einstein-aether theory [38]. 
On the other hand, some important results do not go through from general Einstein-aether 
theory, as they hinge on the exact form of the field equations, and specifically on whether the 
vector is a gradient of a scalar or not. For instance, though qualitatively similar, quantitatively 
the PPN constraints on non-projectable Hofava-Lifshitz gravity cannot simply be inferred from 
those of Einstein-aether theory via a parameter mapping and have been obtain independently 
in Ref. [30]. 



6. Cosmology 

As already mentioned, all versions of Hofava-Lifshitz gravity presented above are written in a 
preferred foliation, and, therefore, they do not constitute diffeomorphism invariant theories, but 
they are invariant under diffeomorphisms that preserve the foliation. This symmetry provides 
enough gauge freedom to choose 

N = l, N* = 0, g ij = a(t)%, (42) 

and bring the ADM line element (5) into the usual Friedmann-Lemaitre-Robertson-Walker 
(FLRW) form 

(, ' J ■ r 2 {d9 2 ■ <hr(J,!o 2 ) 



ds z = -dt 2 + a(ty 



1 — kr 2 



(43) 



under the usual cosmological assumptions of homogeneity and isotropy. Here, a(t) is the scale 
factor and k = +1,0,-1 according to whether the universe is hyper spherical, spatially flat, or 
hyperbolic. 

Clearly, in this setting N has no space dependence. Therefore, its spatial derivatives vanish. 
This will make the actions of the projectable and the non-projectable versions of the theory, eqs. 
(14) and (16) coincide (assuming that all operators allowed by the symmetries of the theory have 
been consistently taken into account). The field equations do not coincide completely, as it does 
make a difference whether N is assumed to be space-independent a priori or whether the latter 
is enforced as a gauge choice. The difference lies on the Hamiltonian constraint and the rest 
of the equations coincide. The former is global in the projectable theory and local in the non- 
projectable theory. So, when it comes to background cosmology, one can study both version 
of the theory at the same time, just taking into account the subtle issue of the Hamiltonian 
constraint. The dynamics of perturbation around the background will, of course, differ in the 
two versions, but we will not discuss perturbations here. 

Cosmology in Hofava-Lifshitz gravity was first studied in Ref. [39, 10] (with detailed balance 
in Ref. [40]). A recent review of cosmology in the projectable version has recently been given in 
Ref. [23], and as explained above, most of the results summarized there are applicable also on the 
non-projectable version. Here we will only give a very brief overview of background cosmology 
and mention in passing some key advantages of Hofava-Lifshitz cosmology. We refer the reader 
to the relevant literature for more details. 

We start by presenting the field equations in an FLRW background. The supermomentum 
constraint is trivially satisfied, so we are left with two equations, corresponding to the first and 
second Friedmann equations. In the projectable version, the Hamiltonian constraint yields 



3A - 1 a 2 V{a) 8vrGArp| 



rfW { ^-Up.- } = 0, (44) 



where 8itGn = M pl 2 , p=—g 1 / 2 5Sm/$N, Sm is the matter action, and 

V(a) = 90 M pl + — + + . (45) 

pi pi 

In the non-projectable case, the Hamiltonian constraint is local and one can do away with the 
integral: 

3A-1 a 2 _ VXo) 8ttGWp 
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The dynamical equations yield in both cases 
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where p = -gV(2/NJg)5S m /8gV. 

Eqs. (44) and (47) and Eqs. (46) and (47) are governing background cosmology in the 
projectable and non-projectable versions respectively. Notice the following: In the non- 
projectable version one can eliminate a from eq. (47) by using eq. (46) to get 

2 a 12a da 3 KH F ' V 1 

Differentiating eq. (46) and subtracting from eq. (48) yields the standard conservation law 

p + 3^-(p + p) = 0. (49) 

This implies that the conservation law for matter is implied by the two (modified) Friedmann 
equations and, thus, imposing it separately is neither needed nor will it constrain the dynamics. 
This is in direct analogy to general relativity. 

Things are quite different in the projectable case, given the global nature of the Hamiltonian 
constraint. It has been argued in Ref. [41] that a global constraint, such as eq. (44), is irrelevant 
for the local patch of the universe inside the Hubble horizon which the FLRW spacetime is 
supposed to approximate. Now, we could, instead of differentiating eq. (44), integrate eq. (47) 
and ignore eq. (44) completely. One then gets 

3A - 1 a* 



V(a) SvrGW / , C(t)\ 
2 a- ^ = ^-( v P+ ^rJ' (50) 

where the functional form of C(t) depends on the conservation law matter satisfies. If we 
eliminate a from eq. (47) by using eq. (50) we get 

3A-1S _ 1 <VM + ^ L ( p+ m +3p Y (51) 



2 a 12a da 3 \ o 3 

Now, the conservation law for matter has become essential to the dynamics via the presence 
of C(t). The terms containing C(t) are the only difference between eqs. (50) and (51), which 
describe the background dynamics in the projectable case, and eqs. (46) and (48), which describe 
the background dynamics in the non-projectable case. 6 



This analysis, leading to the same conclusions, can be performed also before imposing the FLRW ansatz and 
the relevant symmetries, see Ref. [41] . 



Now that we have the field equations in a suitable form for both cases we can proceed to 
discuss phenomenology. Suppose that matter is coupled in such a way that at low energies p 
and p have the usual meaning and satisfy the standard conservation law of eq. (49) (this is a 
strong assumption, see also next section). For the projectable case this implies that C(t) = Co, 
where Co is a constant. The Co related terms then play the role of a dark matter component 
[41]. This dark matter component will not exist in the non-project able case. 

Setting aside the role of matter, the main difference between Hofava-Lifshitz gravity and 
general relativity in background cosmology is the presence of the last two terms in V. The first 
one is what is known as a dark radiation component, as it scales as a~ 4 . The second one is 
what one could call a stiff matter component, as it scales as a~ 6 . Note that the sign of the 
a~ 6 term depends on the sign of the k and that both terms vanish when k = 0. An interesting 
feature related to the presence of these terms is that they can lead to solutions with classical 
cosmological bounces, provided that their coefficients will have appropriate signs [10, 39]. 

One last important property of Hofava-Lifshitz cosmology which is worth mentioning before 
closing this brief overview, is that is seems to lead to a scale invariant spectrum of cosmological 
perturbations, without the need for inflation [42, 43]. This property is related to the anisotropic 
scaling between space and time that the theory exhibits at high energies, which is also the key 
ingredient for power counting renormalizability. 

7. Conclusions and future perspectives 

A brief overview of Hofava-Lifshitz gravity was given and special attention was paid in 
distinguishing between the various different versions. After that, the main focus of this review 
was on the dynamics and the viability of these versions. In summary, the projectable version 
has serious viability problems because it harbours a scalar mode which is either classically or 
quantum mechanically unstable (depending on the position in parameter space) and also exhibits 
strong coupling at low energies (whether or not the Vainshtein mechanism can alleviate the latter 
is not clear yet, see also above). The non-projectable version, on the other hand, seems to have 
similar (and other) problems if detailed balance is imposed. These problems, however, appear to 
be overcome once detailed balanced is abandoned and all operators allowed by the symmetries of 
the theory are consistently taken into account. Alleviating the strong coupling problem requires 
the introduction of a second scale in the theory which should be parametrically smaller than the 
Planck scale, and should be the scale at which higher order operators become important. The 
dynamics of the version exhibiting the extra £7(1) symmetry, which was very recently proposed, 
have not been adequately studied yet. 

There are many open issues regarding Hofava-Lifshitz gravity, the two most important ones 
being the following. Firstly, in the introduction we argued that the theory is power-counting 
renormalizable, as is usually done in the literature. Even though this is a strong indication for 
UV completeness, renormalizability beyond power counting has not been explicitly shown (see 
Ref. [45] for some work in this direction). Additionally, the renormalization group flow of the 
various couplings has not been studied, which implies that, for the time being, one does not 
really know if the theory approaches GR in the IR (A — > 1, r] — > 0) or not. Secondly, the role 
of matter and its coupling to gravity have not been fully clarified yet. The matter action will 
have to include higher order spatial derivatives, which implies that there will be modifications 
in the dispersion relations of matter fields that can lead to serious constraints. Additionally, 
couplings between the matter and the scalar graviton could lead to violations of the equivalence 
principle. Even if initially omitted, such couplings would be generically generated by radiative 
corrections, unless some symmetry prohibits them (e.g. supersymmetry [44]). Hopefully, future 
work will shed some light onto these issues. 
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